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ABSTRACT 

Context. The millisecond pulsar PSR J1903+0327 rotating at 465 Hz has the second highest precisely measured mass 
O . (1.67 Mq) and a weak surface magnetic field (~ 2 x 10 8 G). It is located in the Galactic plane, bound in a highly 

eccentric (e = 0.44) orbit in a binary system with a solar-mass main-sequence star. These observational findings pose 
H ' a challenge for the theory of stellar evolution. 

Aims. Using the intrinsic parameters of PSR J1903+0327 evaluated from radio observations (mass M, rotation period 
P, and magnetic field B deduced from P and P) and a model of spin evolution during the "recycling" phase (spin-up 
by accretion from a low-mass companion lost afterwards) that takes into account the accretion-induced magnetic field 
decay, we aim to calculate the mass of its neutron star progenitor, Mi, at the onset of accretion. In addition, we derive 
■ constraints on the average accretion rate M and the pre-accretion magnetic field B[. We also seek for the imprint of 

r the poorly known equation of state of dense matter at supra-nuclear densities on the spin-up tracks and the progenitor 

neutron star. 

, Methods. Spin-up is modeled by accretion from a thin magnetized disk, using the magnetic-torque disk-pulsar coupling 

O . model proposed by Kluzniak and Rappaport. We adopt an observationally motivated model of the surface magnetic 

field dissipation caused by accretion. We consider three equations of state of dense matter, which are consistent with the 
existence of 2.0 Mq neutron star. Orbital parameters in the accretion disk are obtained using the space-time generated 
by a rotating neutron star within the framework of general relativity. 



Results. Constraints on the progenitor neutron star parameters and the accretion itself are obtained. The minimum 
average accretion rate should be higher than 2 — 8 x 10~ 10 Mq yr _1 , the highest lower bound corresponding to the 
stiffest equation of state. Allowed Bi-dependent values of Mi are within 1.0 — 1.4 Mq, much lower than the oversimplified 
^vq , but widely used B = result, where one gets M; > 1.55 Mq. 

Conclusions. The influence of magnetic field in the "recycling" process is crucial - it leads to a significant decrease in 
■ the spin-up rate and higher accreted masses, in comparison to the B = model. The estimated initial neutron-star 

£f>l ' mass depends on the assumed dense-matter equation of state. We also show that the otherwise necessary relativistic 

, corrections to the Newtonian model of Kluzniak and Rappaport, related to the existence of the marginally-stable 

^s^j ■ circular orbit, can be neglected in the case of PSR J1903+0327. 

\Q . Key words, dense matter - equation of state - stars: neutron - pulsars: general - accretion disks 

o 

1. Introduction through to 2010 January using the 305-m Arecibo radio 

telescope and the 105-m Green Bank Telescope, found that 

, The discovery of PSR J1903+0327, which was the first the pulsar magg ig i. 67±0 .021 M (at the 99.7% confidence 

^ . millisecond radio pulsar found in a binary system with a limit ) and showed strong evidence that effects other than 

»H i 1-03 Mq mam-sequence (MS) companion in an eccentric S p ace _ time curvature, e.g., stellar win ds or tidal fo rces act- 

■ - - ■ ( e = °- 44 ) 95 " da y orblt > poses a challenge to formation ing on Mg companioil; are negligible (|Freire et alll20ll . 
theories of millisecond pulsars (jChampion et al.l . 120081 ). On 

the one hand, its timing parameters - spin frequency / = According to the current theory of neutron-star (NS) 
465 Hz and spin period derivative P — 1.88 x 10~ 20 s s _1 evolution, millisecond radio pulsars (P < 10 ms) originate 
- are nothing unusual for millisecond pulsars. On the other from "radio-dead" pulsars via the accretion- caused spin- 
hand however, the high eccentricity of the orbit and the na- up in low- mass X-ray binaries (LMXB see lAfpar et all 
ture of the companion (main-sequence star), as well as its |1982t iRadhakrishnan fc Srinivasanl 119821 ). This so-called 
location in the vicinity of the Galactic plane make it unique, "recycling" process is believed to facilitate the spin- up 
Pulsar-clock stability and the high eccentricity of the orbit from the initial ~ 0.1 Hz frequency to ~ 500 Hz in 
enabled the pulsar mass measurement to be determined (by ~ 10 9 yrs and is associated with accretion of ~ 0.1 Mq. 
means of the Shapiro delay) with formidable precision. An This idea was corroborated by the detection of mil- 
analysis of the timing data obtained from 2006 December lisecond X-ray pulsations in LMXBs, which were inter- 

preted as the manifestation of r otating and accreting NSs 

Send offprint requests to: bejger@camk.edu.pl ([Wiinands Sz van der Klisl . Il998l) . These millisecond X-ray 
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pulsars are thought to become radio millisecond pulsars af- 
ter the accretion process expires. Radio millisecond pul- 
sars are extremely stable rotators, with P ~ 10~ 20 — 
10~ 19 ss _1 . The surface magnetic fields estimated from 
the timing properties are three - four orders of magni- 
tude weaker than those in normal radio pulsars, for which 
B ~ 10 12 G. This is explained either by the "burying" 
of the original mag netic field under a layer of accreted 
~ 0.1 Mp, material (|Bisnovatvi-Kogan fc Kombera . Il974t 
iTaam fc van den Heuvel Il986l ICumming et all l200lh 
or/and by the Ohmic di ssipation of electric currents in the 
accre tion- heated crust (j Roman 1 119901: iGeppert fc Uroinl . 
Il994ft . The "recycling" in LMXBs is a particularly efficient 
mechanism in dense stellar systems such as globular clus- 
ters and is in accor dance with spe cific statistics of radio 
millisecond pulsars (jLorime E 120081) . Out of a total of 213 
millisecond pulsars, some 100 (i.e., nearly half) arc indeed 
located within binaries, while among the 1850 radio pul- 
sars, only 141 (i.e., 8%) are in binaries. Moreover, out of all 
radio millisecond pulsars, more than half (130) are found 
in 26 Galactic globular clusters! Finally, as much as 40% 
of the 130 globular-cluster radio millisecond pulsars are in 
binaries. 

Efficient spin-up to the kHz frequencies is possible be- 
cause the low-mass companion remains for ~ 10 8 — 10 9 yrs 
in the post-MS phase, overflowing its Roche lobe and feed- 
ing the accretion disk around a NS. Consequently, tidal 
friction within the extended companion has ample time to 
circularize the orbit. Only an additional dynamic perturba- 
tion involving a third star, an occurrence that is not so rare 
within a globular cluster, can make the orbit highly eccen- 
tric or even tear the system apart. This may be the origin 
of eccentric binaries as well as some isolated millisecond 
pulsars, both in globular clusters and the Galactic disk. 

However, PSR J1903+0327 does not fit the above pic- 
ture because neither its orbit is circular, nor its companion 
a post-MS star. It is thus a first specimen of a new group; 
constructing a viable scenario of its formation turns out to 
be a difficult challenge. We briefly and critically review the 
principal scenarios that have been proposed since the dis- 
covery of the PSR J1903+0327 binary in Sect. EJ All but 
one are excluded although the surviving tripl e system sce- 
nario is not free of its own probl ems either (iFreire et all 
120111 : IPortegies Zwart et alll2011|) . Nevertheless, it is clear 
that PSR J1903+0327 was not "recycled" by its current 
MS companion. In the following, we focus on the LMXB 
stage of the pulsar evolution that certainly preceded the 
formation of the presently observed binary. We consider the 
millisecond-pulsar formation model in order to (hopefully) 
deduce the parameters of the pre-accretion progenitor NS 
and to place constraints on the poorly known equation of 
state (EOS) of dense matter. 

We provide a brief introduction to accretion spin- up sce- 
nario and the dissipation of pulsar magnetic field due to ac- 
cretion in Sects. El and [H In Sect. El we describ e the main 
elements of the iKluzniak fe Rappaportl (j2007|) magnetic- 
torque model (hereafter referred to as KR) acting on a NS 
accreting from a thin accretion disk. We report on efforts 
to extend the Newtonian KR model to include the effects of 
the space-time curvature, especially to take into account the 
existence of the marginally-stable circular orbit predicted 
by General Relativity. In our "recycling" simulations, we 
use three EOS of dense matter con s istent with the existence 
of a 2.0 M NS (jDemorest et all l2010j) - these equations 



are briefly described in Sect. El while the results of our sim- 
ulations are presented in Sect. [7] As a final result, we place 
constraints on the average accretion rate during the spin-up 
stage, the initial magnetic field, and the mass of the pro- 
genitor NS. In Sect. El we summarize our main results. An 
Appendix contains a terse overview of the influence of the 
space-time curvature effects, something that was neglected 
in the original KR model, on the spin-up tracks of PSR 
J1903+0327. 

2. Proposed scenarios for the formation of the 
PSR J1903+0327 binary 

Three formation scenarios, denoted below as I- III, were pre- 
vious ly proposed in the discovery paper (jChampion et all 
20081) . and the first scenario was developed further by 
Liu fc Lil (|2009f ). These three scenarios were ruled out by 
new observational data combined with theo retical mod- 
eling and a fourth scenario was advanc ed by IFreire et al.l 
(12011D and iPortegies Zwart et all (|2011h . These scenarios, 
together with their critique, are briefly discussed below: 

/. Rapid rotation at birth 

The pulsar was born spinning rapidly in a core-collapse su- 
pernova, in a binary system in the Galactic disk, with a MS 
companio n, the sup e rnova kick making the orbit strongly 
eccentric (|Liu fc Lil . 12009ft . A spin frequency of 465 Hz 
was reached as a consequence of accretion from the fall- 
back disk, which requires less than 10 3 — 10 4 yrs before 
the accretion (or rather hyperaccretion - the mean rate is 
~ 10~ 4 M Q yr _1 ) stops and the pulsar becomes an ac- 
tiv e radio emitter. This sc enario was considered unlikely 
by IChampion et al.l (|2008f) , because of the very low sur- 
face magnetic field (B ~ 10 8 G) implied by the measured 
value of P. iLiu fc Lil 1)20091) proposed that the low magnetic 
field value was caused by its dissipation associated with ac- 
cretion. However, they used a model of a magnetic field 
dissipation developed originally for slow accretion in an 
LMXB, which is not valid for a hyperaccreting newly born 
NS. Additional arguments, in partic ular observational ones, 
against this scenario can be fo und in lFreire et al.l (1201 lh and 
IPortegies Zwart et al.l (|2011l ). which is therefore unaccept- 
able. 

//. Hierarchical triple system 

The pulsar was recycled in a LMXB and is currently a mem- 
ber of a hierarchical triple system containing a white dwarf 
(inner binary with the pulsar) and a MS star (outer binary) . 
The high eccentri city of the pu lsar orbit is generated by the 
Kozai resonance (jKozail Il962l ) between the inner and outer 
binaries. The scenario is ruled out by the precise measure- 
ment of the optical spectra of the MS companion of the 
pulsar permitting the calculation of the companion radial 
velocity, which turned out to be consistent with predictions 
based on the orbital parameters determined from the pulsar 
timing. The optical data lead to an independent estimate 
of the mass ratio in the binary (jFreire et all 120111 ) . Finally, 
there was no time dependence of the eccentricity e in the 
TEMP02 timing analysis, while e growth is expected ow- 
i ng to the Kozai resonance in the tr iple hierarchical system 
dGopakumar. Bagchi fc Ravi I2009D . The scenario should 
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ther efore be ruled out as inconsistent with observat ions (see 
also iFreire et alJl201lt iPortegies Zwart et al]l20ll 



///. Pulsar spun up in a LMXB and then ejected into the 
Galactic disk together with a newly captured MS companion 

The pulsar was spun up in an LMXB within a globular clus- 
ter and then exchanged the evolved companion in a process 
of interaction with a MS star in the cluster core, its binary 
orbit becoming very eccentric. The kick-off in the three- 
body interaction was sufficient to eject the final NS-MS bi- 
nary out of the globular cluster and into the Galactic disk. 
However, measurements of the proper three-dimensional ve- 
locity of_tLe__bJ^ry_allowed to track its position back in 
time (|Freire et all 20111) . The binary was found to always 
be within 270 pc from the Galactic disk and farther than 
3 kpc from the Galactic center. Therefore, an exchange in- 
teraction (where an evolved companion in a LMXB is ex- 
changed for a MS star), which is probable in a dense stel- 
lar environment such as a globular cluster or the Galactic 
center, is so unlikely that i t leads to the rejection of the 
scenario (|Freire et all 1201 it ). 



IV. Pulsar in a triple (tertiary) system with two MS stars of 
different masses, the pulsar being spun-up by accretion from 
the evolved, more massive companion that was afterwards ex- 
pelled from the tertiary to make it a binary 



This scenario was proposed by IFreire et al.l (|201lD and 
IPortegies Zwart et al.l (|2011l) . A starting point is a tertiary 
composed of a massive star and two MS stars of differ- 
ent masses. The massive star collapses, giving birth to a 
pulsar, the tertiary still being bound. The initially more 
massive MS star (MSI) evolves forming a LMXB with the 
pulsar and spinning it up by accretion. After a substan- 
tial mass loss, MSI is removed from the system owing to 
the ablation/accretion induced by the interaction with the 
millisecond pulsar or ejected fr om the system by means 
of the three-body interactions (jFreire et all 1 2 1 lh . This 
leaves a system as the observed one, with a 1.03 M Q com- 
panion still at the MS stage, in an eccentric orbit. The 
ablation/accretion removal of the companion was previ- 
ously proposed as a mechanism producing isolated recy- 
cled pulsars in the Galactic disk. Ablation by the NS radi- 
ation indeed operates in the PSR B1957+20 binary (Black 
Widow), whose companion mass was reduc ed to 0.025Af Q 
(jFruchter et all 119881: iKluzniak et~aH 119881) . However, de- 
tails of the ablation process are rather uncertain and the 
timescale needed to completely vaporize the compan i on ca n 
be longer than the Hubble time (|Levinson fc Eichleit 1 1 9 9 J) . 
Alternatively, the expulsion of the lower-mass companion 
owing to the chaotic character of a three-body interaction in 
the triple syste m is proposed. We n ote, however that refer- 
ences quoted in IFreire et al.l (|201ll) do not see m to provide 
an accurate explanation of PSR J1903+0327: iHutl (1 1984) 
assum es equal masses for the three bodies, and iPhillipsi 
(|l993h studies the reflection of the pulsar radio beam from 
the swarm of asteroids orbiting the pulsar and derives an 
upper bound of ~ 10~ 4 M to the mass of the asteroid-like 
material orbiting the Vela pul sar within 1 AU. In our opin- 
ion, the scenario presented bv IFreire et al.l (|201lf ) is not yet 
properly supported by relia ble quantitative estimates. In 
IPortegies Zwart et al.l (|201ll ). both the evolution and the 



dynamics of the NS+MS1+MS2 system are studied in de- 
tail. Several competing channels of evolution for the ter- 
tiary are compared using advanced numerical simulations. 
The available parameter space leading to the observed PSR 
J1903+0327+MS2 binary is estimated and the birth-rate of 
these binaries in the G a lactic disk is found to be acceptable. 
Portcgics Zwart et al.l ((201 lh highlighted the need to (ob- 
servationally) find a ,, missing link" in the scenario: a wide 
LMXB orbited by a tertiary low-mass MS star. It is clear 
that the discovery of such a system would strongly fortify 
the triple-star scenario. 

Four scenarios, reviewed critically above, have been pro- 
posed until now. Only one of them cannot be immedi- 
ately ruled out by observations, but even this one lacks 
a solid quantitative basis. However, it is clear that PSR 
J1903+0327 has been "recycled" in order to explain its 
present spin period, low surface magnetic field, and mass, 
which is much higher than the average 1.4 Mq. From now 
on, we focus exclusively on the spin-up phase of the pulsar 
evolution, to deduce its pre-accretion parameters. 

3. Spin-up by accretion in LMXBs 

We consider a binary consisting of a young radio pulsar (of 
a "canonical" surface magnetic field B\ ~ 10 12 G, and a 
spin period of a fraction of a second) and a low-mass MS 
companion. In a few million years, it spins down by means 
of magnetic dipolc braking down to the period of a few sec- 
onds, without suffering any significant magnetic field decay. 
Consequently, it crosses the radio pulsar death line and dis- 
appears as a pulsar. On a much longer timescale of a billion 
years, the low-mass companion enters the red giant phase 
and fills its Roche lobe. This entails the mass transfer onto 
the NS via an accretion disk. The binary system becomes a 
LMXB, remaining in this stage for 10 s — 10 9 yrs. Accretion 
of the plasma onto a NS increases its mass, while accelerat- 
ing its rotation, as well as inducing the decay (dissipation) 
of its surface magnetic field. 

The accretion rate is expressed in terms of the baryon 
mass Mb, since it is the parameter that can be uniquely 
determined for the binary system; indeed Mb is related to 
the star's total baryon number iVb by Mb = N^mo, where 
mo is the mass per nuclcon of the 56 Fe atom. The increase 
in Mb is straightforward and proportional to the accretion 
rate - at time t (measured by a distant observer) it is de- 
noted as Mb(t). Assuming that accretion starts at t;, the 
integrated (total) baryon mass increase is 



AMb(t) 



M h {t')dt'. 



(1) 



Since the detailed history of accretion is unknown, we treat 
M b as a constant, such that AM h (t) w M b (i— While M b 
is a very important global stellar parameter, the quantity 
that is actually measured is the gravitational mass M. We 
define its increase as 



AM (t) = J M{t')dt'. 



(2) 



It is found that AM depends on AMb according to the re- 
lation betw een dM, dMb, and the total angular momentum 
J (see e.g.. lFriedman et aT1ll988l ) 



dM = ildJ + u'dM b 



(3) 
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where u , the time component of fluid 4- velocity as seen by 
a distant observer, has the meaning of chemical potential 
(divided by m c 2 ). For AM(i), one then obtains 



AM(t.) = J M b («* + i tot n)dt' J 



(4) 



using the total angular momentum J evolution equation 



dJ 
dM b 



l tot (M h ,f,B), 



(5) 



where l to t includes the angular momentum transferred to 
the star by the infalling matter and the influence of the 
magnetic torque (for details, sec Sect. [5]). The values of 
the gravitational mass M and spin frequency / = 1/P for 
a given Mb and J are calculated f or stationary rigidly ro - 
tating two-dimensional NS models (|Bonazzola et al.l . ll993|) . 
with the rot star code implementation from the numerical 
relativity library LORENEfl 



4. Decay of neutron-star magnetic field in LMXBs 

Millisecond pulsars have typically a low surface po lar mag- 
netic field of strength B p ~ 10 8 — 10 9 G (see e.g.. lLorimerl 
I2008T ) . but observations do not provide any evidence of 
the Bp decay during the radio-pulsar phase. However, a 
substantial B p decay (of some four orders of magnitude) 
is thought to occur during the accretion "recycling" in a 
LMXB, leading to the format ion of a millise cond pulsar 
(|Taam fc van den Heuvellll986l for a review see lColpi et alJ 
|2001[ ). We relate B p to measured values of P and P 

, , \ 1/2 

by a standard formula, B p = 3.2 x 10 19 I PP/s) G 

(|Manchester fc Tavlorl Il977t ). This assumes a stellar mo- 
ment of inertia / = 10 45 g • cm 2 and radius R = 10 km, 
and enables one to replace the measured value of P by B p . 



4.1. Decay of magnetic field - theory 

Theoretical modeling of accretion-induced decay of B p 
turns out to be a challenging task. The original idea 
that B p decays in close binaries because it is 'buried' 
('screened') by the accreted matt e r was p r opose d by 
iBisnovatvi-Kogan fc Komberd (| 19741 ). [Roman i (199Q|) sug- 
gested that a combination of crustal heating due to accre- 
tion, accelerating Ohmic dissipation, and advection of mag- 
netic field lines from the poles towards the equator could 
explain the B p decay. Some other authors focused on the 
realistic modeling of the acceleration of the Ohmic dissi- 
pation of the crustal B caused by the heating induced b y 
accretion (jGeppert fc Urpinl[l99l lUrpin fc Geppertlll995l) . 
The decay of B p due to the di amagnetic screenin g by 
the ac c reted layer was studied in IZhang et alJ (|l994f ) and 
I Zhang! (|l998f ) . Dependence of the screening of crustal mag- 
netic field o n the accretion rate w as studied in more de- 
tail later bv lCumming et all (|200ll ). These authors consid- 
ered in detail the interplay between the advection of mag- 
netic field and its Ohmic diffusion, and used realistic micro- 
scopic models of the outer 'ocean' (molten crust) and the 
crust. Magneto-hydrodynamical simulation of accretion- 
burial of -Bp have also been performed by several authors 



fe.g. IPavne fc MelatosI I2004L 120071 see also the study of 
IWette et alJ 2010l and references therein). This very brief 
and incomplete (especially as far as the references are 
concerned) review illustrates the theoretical effort to ex- 
plain the decay of B p during the accretion-driven recy- 
cling of the millisecond pulsars. Unfortunately, no reliable 
and robust scenario based on realistic microphysics com- 
bined with complete magneto-hydrodynamical treatment, 
and consistent with astronomical observations, is available 
today. Therefore, in what follows we limit our description 
of the Bp decay to simple phenomenological models, based 
to some extent on observations of three populations: old ra- 
dio pulsars, accreting binary neutron stars, and millisecond 
pulsars. Our basic assumption (which is widely accepted) 
is that these populations are linked by the recycling pro- 
cess, which in turn is driven by accretion onto a magnetized 
neutron star in a close binary. 

4.2. Decay of magnetic field - phenomenology (with some 
observational basis) 



iTaam fc van den Heuvell (| 1986T) analyzed a set of LMXBs 
of different ages and therefore different amounts of ac- 
creted mass. They suggested that there is a possible in- 
verse correlation between B p and the (estimated) total 
amount of accreted material. Their conclusion was con- 
firmed in a later study by Ivan den Heuvel fc Bitzarakil 
(|1995l ). IShibazaki et al.l (|1989l ) presented more detailed ar- 
guments based on a subset of LMXBs that enabled this 
inverse correlation to be quantified as 



Bp = B p (AM h ) = A/(l + AM h /m B ), 



(6) 



1 http://www.lorene.obspm.fr 



where B\ is the initial (pre-accretion) magnetic field, 
Bp(AMb) is the magnetic field after AMb = M^t is ac- 
creted by the neutron star, and me is a constant setting 
the scale of dissipation of B D with increasing AMb. We 
note that Shibaza ki et alJ ([1989) do not advocate any phys- 
ical model leading to Eq. ([6]). They report astrophysically 
interesting bounds on me, resulting from the application 
of Eq. © to some LMXBs, X-ray pulsars, and millisecond 
pulsars. In particular, for me ^ 10~ 3 Mq, the evolution 
at an X-ray pulsar stage proceeds along the equilibrium 
spin-up line, with accret ion (advection) torqu e balanced 
by the magnetic torq ue (jGhosh fc Lambl fl979l Fig. 1 of 
IShibazaki et al.lll989h - rapid millisecond rotation cannot 
be obtained in this way. If, however, me 10 -4 Mq, 
then the evolutionary tracks in the log-B p — logP plane 
are consistent with measured pairs of B p and P for bi- 
nary and isolated millisecond radio pulsars. Finally, tak- 
ing tub ~ 10~ 4 Mq allows a reasonable representation 
of the inverse correlation betw een B p and AMb noted by 
ITaam fc van den Heuvel (|l986l 

One has to be aware that Eq. ^ is based on a 
limited and uncertain set of data referring to LMXBs. 
It is clear that Eq. ^ is too simplistic to describe 
magnetic field decay for all kinds of accreting neu- 
tron stars. It is natural to expect that the decay of 
Bp depends not only on AMb, but also on M, which 
i s decis i ve in the heating of t he neutron-star interior 
(IWiiersl. Il997t lUrpin fc Geppertl . fl99l lUrpin et all . 119981: 
ICumming et al.l . l200lh . Specifically. IWiiersI (|l997l) demon- 
strated that the decay law B p ex 1/AMb is inconsistent 
with a broader set of available data for accreting neutron 
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stars in both X-ray binaries and recycled millisecond pul- 
sars. 

In spite of the limitations and uncertainties discussed 
above, we adopt Eq. ([5]) as our baseline description of the 
Bp decay in LMXBs and compare it with other proposed 
phenomenological formulae for the accretion-induced B p 
decay to assess the model-independent features of our basic 
results. 

For example. IWiie rs (1997]) discussed the square depen- 
dence on AMb /m-Q in the denominator of Eq. ^ as well: 



B p = S p (AM b ) = Bi/(1 + AM h /m B ) 



(7) 



In this case, one has to use a lower value of m-Q ~ 1CP 3 -j- 
10 -2 Mq, which allows to form a millisecond pulsar before 
the field decays to very low values. 

An exp onential decay o f B n w as employed bv lKiel et al.l 
(|2008t ) and lOsIowski et all (|201l[) to be 



B p = (S; - B min ) exp(-AM b /m B ) + B u 



(8) 



where B min = 10 8 G is the assumed minimal residual mag- 
netic field, to reproduce the observed P — P distribution by 
mean s of the population synthesis studies (jOstowski et al.l . 
1201 lh . We discuss briefly how the final results depend on 
the assumed magnetic field decay model in Sect. 17.31 

5. Spin-up by disk accretion onto a magnetized 
neutron star 

We assume that the evolution of an accreting NS can be 
represented as a sequence of stationary rotating configura- 
tions of increasing baryon mass. We use the KR formalism 
to determine the circular orbit ro, from where the accre- 
tion effectively takes place - the radial inner boundary of 
the Keplerian accretion disk in which the viscous torque is 
non-vanishing. Following KR, we define the corotation ra- 

1/3 

dius r c ee (GM/uj 2 ) ' (i.e., the radial distance at which 
the Keplerian orbital angular frequency Qk is equal to 
the rotation frequency of a central star, uj s = 2irf), the 
magnetosphcric radius r m ee (GAf) -1 ^ 7 



M" 2 /V /7 (where 



B p RJ 



approximates the stellar magnetic dipole mo- 



ment), the magnetic-to-corotating radius ratio £ = r m /r c , 
and the so-called "fastness" parameter u = uj s /flK{ro) = 

(ro/r c ) 3 / 2 . In addition, we include relativistic effects - in- 
evitable in the space-time generated by a rotating NS - 
especially the existence of the marginally stable orbit r ms 
and denote the Schwarzschild gravitational radius as r s = 
2GM/c 2 , and the corotation radius in r s units, /3 = r c /r s . 

We then solve the equation that corresponds to the van- 
ishing of the viscous torque to determine the disk inner 
boundary, namely 



dr 



Mr 4 



(l-<0, 



(9) 



where I is the specific angular momentum of a particle cal- 
culated in a relativistic way. To avoid the arduous calcula- 
tions of the NS space-time that would be needed to obtain 
the specific angular momentum I at a given disk radius r 
in every time-step o f the evolution, we exploit the result of 
iBeieer et al.l (|2010T ) . who provide an approximation of the 
value of l(r), defined as 



l(r) 



y/l - V 2 / C 2 



(10) 



with an approximate value of the particle orbital velocity v 



I GM 2GJ 



y/l - 2GM/(rc 2 ) \ V r 3 r 3 c 2 



(11) 



This approximate value of /(r) deviates by less than one 
per cent from the true value for spin frequencies and masses 
similar to PSR J1 903+0327 measurements (for details sec 
IBeieer et al.ll2010l) . Eqs. (pi fTTj) yield an algebraic equation 
for ro of 



1 



/ms(ro) = 




(1-w), (12) 



with a dimensionless function / ms defined as 

f (r) - — — 
s lr dr 

which for l(r) given by Eqs. (fT0)l and ((TT1) equals 

1 - a/f 3 / 2 



(13) 



/ms(fo) 



(l-V 2 /c 2 ) 3 / 2 y/l-l/f 

3a 



r — 2 v 
x I 2 



1 



F 3/2 



(14) 



where a = Jc/(y/2GM 2 ) and f = r /r s = /3lu 2 / 3 . 

In general, Eqs. (|12p and (fTfj) allow us to calculate ro 
in the "recycling" process of even very rapidly rotating and 
massive millisecond pulsars, when the existence of a rel- 
ativistic marginally-stable orbit cannot be neglected i.e., 
when the KR condition of ?*o <C r ms is no longer valid. The 
value of r ms corresponds to the solution of 



dl 
dr 



0. 



(15) 



which is determined by assuming that / ms (r ms ) = 0. The 
validity of employing this general, refined approach in the 
case of PSR J1903+0327 is studied in the Appendix. 

To calculate the increase in the total stellar angular mo- 
mentum J, we take into account the transfer of specific an- 
gular momentum Iq ee /(ro) and the KR prescription for the 
magnetic torque 



dJ 
dAih 



= hot = K r o) - 




9r 3 M b 



(16) 



6. Equations of state 

The EOS of dense cores of NSs remains poorly constrained. 
This is due to, on the one hand, a lack of knowledge of 
strong interactions in dense matter, and on the other hand, 
deficiencies in the available many-body theories of dense 
matter. This uncertainty has been reflected as a rather 
broad scatter in the theoretically derived and EOS de- 
pendent maximum allowable masses for NSs, M max (EOS) 
(see, e.g. lHaensel et al.l [20071) . Fortunately, the measure- 
ment of the mas s of P SR J1614-2230, of 1.97 ± 0.04 M Q 
(jDemorest et al.l . I2010T ) . introduces a rather strong con- 
straint of M max > 2.0 M.q. This means that the true 
EOS is rather stiff. To illustrate a remaining uncertainty 
in the stiffness, we considered three different models of the 
EOS. In all cases, the simplest composition of matter was 
assumed - neutrons, protons, electrons, and muons in /?- 
equilibrium (npefi): 



6 



Bejger et al: Neutron star progenitor of PSR J1903+0327 



DH This model of iDouchin fe Haensell (|200lt ) is non- 
relativistic and its energy density functional is based 
on the SLy4 effective nuclear interaction. The model de- 
scribes in a unified way both the dense liquid core of NS 
and its crust, yielding M max = 2.05 M Q and a circumfer- 
ential radius at the maxim um mass of R m = 10.0 km. 

APR = Al8+5v+\JlX* from lAkmal et al.l (jl998t ) is a non- 
relativistic model with some relativistic corrections. It 
consists of a two-nuclcon Argonne potential A18 with 
relativistic boost corrections Sv and an adjusted three- 
nucleon Urbana UIX* potential. A variational solution 
of the many-body problem yields M max = 2.21 M Q 
and a circumferential radius at the maximum mass of 
R Mmax = 10.0 km. 

BM = TM16S0 (with some minor changes) was drawn 
from a set of relativistic models of iBednarek fc Mahkal 
(|2009t ). It consists of a Lorentz-covariant effective non- 
linear Lagrangian including up to quartic terms in me- 
son fields, based on chiral-symmetry breaking expan- 
sions. The EOS was calculated in the mean field approx- 
imation, yielding M max = 2.11 M and a circumferen- 
tial radius at the maximum mass of i?M max = 11.95 km. 




R [km] 



7. Results 

To illustrate the laws and relations governing the ,, recy- 
cling" process, we begin by constraining the final spin fre- 
quency, / = 465 Hz. We then calculate sets of the evolu- 
tionary tracks labeled with the NS initial parameters - A/;, 
Pi, and Bi - covering a broad range of possible values. If 
not stated otherwise, the figures relate to results obtained 
using the DH EOS. In what follows, we denote B p by B 
and M b by M. 

7.1. Constraining the parameter space: final frequency 

f = 465 Hz 

Fig. Q] presents the mass-radius relation for accreting NSs 
for different values of accretion rate M and initial magnetic 
field Bi. The initial mass is fixed at 1.4 M©. The upper ends 
of each curve correspond to the final frequency, 465 Hz. 
For Bi = 10 12 G, the lowest accretion rate considered is 
3 x 10~ n M©/yr because for lower rates the configurations 
enter the axisymmetric-perturbation instability region al- 
ready for frequencies lower than 465 Hz. The same configu- 
rations are shown in Fig. [2J where the rotational frequency 
as a function of mass is presented for different accretion 
rates. For some evolutionary tracks in Fig. [2J the initial 
frequency is set to fa = 50 Hz, the value corresponding to 
a ty pical frequency expected for new ly-born pulsars (Table 
7 of iFaucher-Giguere fc Kaspl 1200(1 . These curves are in- 
distinguishable from those corresponding to initially non- 
rotating configurations, as the braking timescale is by many 
orders of magnitude shorter than the spin-up time. The dif- 
ference is visible on a logarithmic scale in Fig. [3J where we 
present the frequency evolution of an accreting star for two 
initial configurations: non-rotating NS and a configuration 
rotating initially with frequency 50 Hz. The accretion rate 
here is fixed to til = 10~ 9 M /yr. For B\ ~ 10 12 G, the 
spin-down predominates for ~ 10 3 yrs, and afterwards the 
accretion spins-up the star. The amount of accreted mate- 
rial depends very sensitively on the strength of the initial 
magnetic field; correspondingly, for Bi = 10 11 G the spin- 



Fig. 1. Mass-radius relation for accreting stars with differ- 
ent initial magnetic field Bi and accretion rate M . Solid 
black curve denotes static configurations. Evolutionary 
tracks (arrows mark the direction of evolution) correspond 
to the following initial parameters (Bi [G], M [M Q /yr]) - 
a: (10 12 , 3 x 10" 11 ), b: (10 12 , 10" 10 ), c: (10 12 , 10" 9 ), d: 
(10", 10- 11 ), and e: (10 11 , 10~ 9 ). Tracks c and d coincide, 
as explained in Sect. 17. ll (color online). 
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Fig. 2. Spin frequency evolution during accretion, for ini- 
tial frequencies zero and 50 Hz, as a function of stellar mass. 
Three different cases of accretion rates (in M Q /yr) for ini- 
tial magnetic field B t = 10 12 G (solid line), 10" G (dashed 
line) are shown (color online). 
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Fig. 3. Spin frequency evolution for initial spin frequen- 
cies of zero and 50 Hz and for an initial magnetic field of 
Bi = 10 12 G and 10 11 G as a function of the accreted mass 
(upper axis) or time (lower axis) calculated for a constant 
accretion rate M = 10~ 9 M Q /yr, on a logarithmic scale 
(color online). 



down timescale is two orders of magnitude longer (10 5 yrs) 
than for B{ = 10 12 G. This is a direct consequence of 
the quadratic dependence of the magnetic torque on the 
magnetic field B, Eq. (|16[) . The minimum value of rota- 
tion frequency corresponds to the exact balancing of an- 
gular momentum ^o(ro) at the accretion disk edge, To, by 
the magnetic torque. In Figs. 0] and we show the mag- 
netic field vs. mass dependence in such a case (the error 
bar corresponds to the 3 — a measurement of the mass, 
M = 1.667 ±0.021 M ). 

We note that the problem under consideration is degen- 
erate with respect to the constant 



q = B 2 /M, 



(17) 



(sec Fig.0]), as a consequence of Eqs. (fT2"l) and (|TB)) depend- 
ing on the quantity q, and not B and M separately. 

To obtain the results for some M%, the magnetic field 
B corresponding to M\ should be multiplied by the factor 



M2/M1. This relation allows us to determine the lower 
bound to the accretion rate for given observed values of B 
and M or, assuming an accretion rate, to determine the 
maximum value of the final magnetic field e.g., if M = 
M w = 10~ 10 M /yr (red curve on Fig. HJ, 5 max = 1.16 x 
10 s G for M = 1.67 Mq. From the observations, we infer 
that Bf ~ 2 x 10 s G; since the condition £? max > Bf should 
be fulfilled, we use the scaling law 



B m a X (M) = B 

max (M10) 



M 



1/2 



M- 



(18) 



Fig. 4. Final magnetic field Bf versus gravitational mass 
M for configurations rotating at / = 465 Hz. Colors cor- 
respond to different accretion rates. For a given accretion 
rate (in MQ/yr), different curves are defined by a given 
initial configuration (i.e., central density or initial mass). 
Along each curve, the initial magnetic field increases up- 
wards. The error bar corresponds to the uncertainty (3a) 
in the mass measurement (color online). 



to obtain a lower limit to the accretion rate 



M > 



B f 



Mio = 3 x 10~ 10 M Q /yr. 



(19) 



In Fig. [5J we present the relation between the magnetic 
field and the gravitational stellar mass for the configuration 
spun up to 465 Hz, for three different EOSs described in 
Sect. El The lower limit to M assessed above is marginally 
consistent with the DH EOS for the measured values of 
M = 1.67 M m and B = 2 x 10 8 G. 



7.2. Constraining the parameter space even further, f = 465 
Hz and M= 1.67 M Q 

Lines presented in Fig. correspond to different initial pa- 
rameters (magnetic field Bi and mass Mj) that lead to 
configurations of M = 1.67 M at the spin frequency 
/ = 465 Hz. We employ two fixed accretion rates of 
M = 10~ 10 M Q /yr and M = 10~ 9 M Q /yr. The scaling 



10 



relation Bf(Mi) = B f (M 2 ) ■ ^M 1 /M 2 holds and allows us 

to set limits on M for a given value of the magnetic field. 

The differences between the DH and APR EOSs are 
small in comparison to the difference of them both with a 
much stiffer BM EOS. The value of the stellar moment of 
inertia is responsible for such a discrepancy: I(M) for the 
BM EOS is about 25% higher than the corresponding values 
for the DH and APR EOS (for the same M). Consequently, 
to obtain the same final frequency one needs a larger J and 
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M, [M ] 

Fig. 5. Final magnetic field Bf vs the gravitational mass for 
a star rotating at / = 465 Hz. Different colors correspond to 
different EOSs (solid red - DH, dashed blue - APR, dotted 
green - BM). The accretion rate is M = 3 x 1CT 10 M Q /yr. 
The error bar reflects the uncertainty (3er) in the mass mea- 
surement (color online). 
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Fig. 6. Final magnetic field, Bf, versus initial gravitational 
mass, for a star with final parameters M = 1.67 M Q and the 
frequency / = 465 Hz, for APR, DH, and BM EOSs. The 
accretion rate is M = 10~ 10 M Q /yr (solid curves) and M = 
10 -9 M Q /yr (dashed curves). Scaling relation Bf(Mi) = 

Bf{M 2 ) ■ \J Mi /M2 is fulfilled, i.e., dashed curves match 
the solid ones multiplied by ylO (color online). 




1.2 1.4 1.6 



M [MJ 

Fig. 7. Spin-up tracks of the accreting NSs leading to the 
final configuration rotating at / = 465 Hz and with gravita- 
tional mass M = 1.67 M Q and magnetic field B = 2 x 10 s G. 
Curves arc labeled by the average accretion rate (in Mq /yr 
and initial value of the magnetic field (in G) . For compari- 
son the spin- up for B = (dashed line) , via accretion from 
marginally stable orbit, is shown (color online). 

a larger amount of accreted mass, which leads to the lower 
value of Bf. 

7.3. From a progenitor NS to PSR J1903+0327: 

M = 1.67 M , / = 465 Hz, and B = 2 x 10 8 G 

One can obtain a stringent limit on the parameters of the 
progenitor NS by fixing the final value of the magnetic field. 
Some tracks considered in this subsection are shown in Fig. 
labeled by the initial value of the magnetic field Bi and 
the accretion rate M . The final magnetic field is fixed to 
B = 2 x 10 8 G. For comparison, we present the case in 
which we neglect the effect of the magnetic field entirely: 
the accretion takes place from the marginally-stable orbit 
calculate d exactly in accord ance with General Relativity 
(see, e.g.. lZdunik et al.ll2002t) . 

The spin- up in the B — case is very efficient, one needs 
~ 0.1 Mq only to reach the observed frequency of 465 Hz. 
Neglecting a magnetic field therefore leads to a rather high 
progenitor-NS mass, 1.58 Mq, whereas reasonable values of 
the progenitor B and M enable Mf to be reached starting 
from a moderate value of Mj ^1.3 Mq. The amount of 
accreted matter in the B = case can be larger if we al- 
low for the reduction in the ang ular momentum tra nsport 
efficiency (parameter x\ <C 1 in IZdunik et al.|[2002[ ). here 
however we assume a 100% efficiency (xi = 1), in accor- 
dance with recent numer ical calculations dBeckwith et all 
12001 iShafeeetall 12001 . 

The accretion rate as a function of Mi (assuming final 
Bf = 2 x 10 8 G) is presented in Fig. [SI The lower limit 
to the accretion rate depends on the EOS and ranges from 
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Fig. 8. Average accretion rate vs initial mass needed to Fig. 9. Accretion rate versus time of accretion needed to 
evolve a NS by disk accretion to the frequency / = 465 Hz, reach observable configuration (color online), 
mass M = 1.67 M , and magnetic field B = 2 x 10 8 G for 
the three EOSs under consideration (color online). 

8. Conclusions 



2.5 x 1CT 10 M Q /yr (DH, APR) to 8.5 x lCT 10 M Q /yr (BM). 
For a given EOS, the required M depends weakly on the 
accreted mass needed to spin up the star to 465 Hz, pro- 
vided that AM > 0.2 M . For AM < 0.1 M Q , the required 
accretion rate increases dramatically with decreasing AM. 
Using AMb = Mbt, we obtain the time needed to spin up 
the star to its presently observed frequency. The result is 
presented in Fig. [5] 

We also performed calculations for other B(AM) de- 
pendences (Eqs. [7] and [5] in Sect. @|. The results arc qual- 
itatively similar in the case of an inverse quadratic depen- 
dence (Eq. [7]); the value of M at the "plateau" region in 
Figs. [5] and H is slightly higher, 3.6 x 10~ 10 M /yr, than 
3 x 10" 10 M Q /yr, while the Eq. © is adopted. For the ex- 
ponential field decay (Eq. [5]) however, the resulting M is 
8 x 10~ 10 M /yr (f or the mass dec ay scale tub = 0.05 M , 
which according to lOslowski et all gives the closest agree- 
ment between their population synthesis model predictions 
and the observed P — P data for Galactic pulsars). This 
discrepancy is due to the different behavior of B p (AM) 
close to the final B p value, where Eq. © gives significantly 
larger values of B p than Eq. ^j. For an efficient spin- up, 
the "magnetic term" proportional to B p /M in l to t should 
be small (Eq. ITB|) , hence the relatively large value of B p has 
to be compensated for by a higher value of M - for exam- 
ple, assuming B { = 10 12 G, Eq. © yields A M = 0.5 M w in 
1.7x 1 9 yrs, while the exponential decay of lOslowski et al.l 
(|2011h gives 0.46 M in 6 x 10 8 yrs, that is, the accretion 
of almost the same amount of matter but somewhat more 
rapidly. One should however treat this kind of empirical 
formulae with caution, since in the case of parameters ob- 
tained by population synthesis methods they may be un- 
controllably affected by other assumptions. 



Our simulations have allowed to estimate the intrinsic pa- 
rameters of the progenitor NS required to reach, in the 
process of recycling, the measured parameters of PSR 
J1903+0327. To some extent, as we have shown, the pro- 
genitor NS parameters depend on the EOS of dense mat- 
ter. This was studied usin g three EOSs consisten t with the 
measurement of 2M NS (|Demorest et all l2010l) . 

We have found that the mean accretion rate M during 
recycling should be larger than (2.5 -=- 8.5) x 10~ 10 M /yr, 
the highest lower bound being obtained for the stiffest EOS. 

For each EOS, the required mean accretion rate is ap- 
proximately constant for a broad range of initial masses, 
1 M -=-1.4 M . Therefore, depending on the initial 
magnetic field Bi we can reproduce parameters of PSR 
J1903+0327, or more generally, observed millisecond pul- 
sars, for a specific range of initial masses. In other words, 
the present parameters of a recycled millisecond pulsar do 
not allow us to determine its initial mass, in contrast to 
the case of recycling by accretion with B = 0. Simulations 
that n eglect the magnet ic field (following the seminal pa- 
per of iCook et all Il99l give a rather high value of the 
lower bound to the progenitor NS mass of Mj > 1.55 M 
(1.58 M Q for DH and APR EOS, and 1.55 M for the 
BM EOS). This is a direct consequence of the finding that 
the magnetic field significantly decreases the spin-up rate 
(effectively decreasing the efficiency of angular momentum 
transfer onto the star, as in to the B — and case). 
Accounting for the magnetic field effect (magnetic torque) 
is therefore crucial for the understanding of the observ- 
able pulsar population properties, especially in view of the 
proposal that gravitational wave emission is a major dis- 
sipat i ve agent that prevents efficient spin-up (jArras et all 
[200l IWatts et all [20081 iWatts fc Krisrmarll2CI0^ In other 
words, for evolutionary reasons there may not be sufficient 
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Fig. 10. Upper panel: spin-up tracks of accreting NSs lead- 
ing to presently observed PSR J1903+0327 parameters 
on radius-spin frequency plane. The relevant characteris- 
tic radii are: stellar radius (solid curves), the radius of 
a marginally stable orbit (dotted curves), and the radius 
of the inner boundary of accretion disk (dashed curves). 
Colors correspond to the specific values of accretion rate 
and the initial value of the magnetic field. Lower panel: ef- 
fective radial potential V (r, l(r)) in the energy (mc 2 ) units 
for an exact, approximate, and Schwarzschild solution; the 
dot marks the final tq. See the text for more details (color 
online). 



time and/or a sufficient amount of matter to accrete, to 
form a rapidly-spinning pulsar with a magnetic field suffi- 
ciently strong to produce a detectable radio beam. 

The framework presented here is suitable for testing 
the global properties of the Galactic pulsar population, as 
well as for studying other millisecond pulsars with pre- 
cisely measured masses, e.g., massive millisecond pulsars 
recycled in intermediate-mass X-ray binaries, where the re- 
cycling proce s s is m uch shorter than the one studied here 
(jTauris et all 1201 lh . Our analysis can also be easily ex- 
tended to take into account time-varying accretion rates 
and more detailed multi-parameter magnetic field decay 
prescriptions - this paper is the first of a series devoted 
to these studies. 
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Appendix: The importance of relativistic effects in 
the "recycling" process of PSR J1903+0327 



The original framework of Kluzniak & Rappaport ( 2007) is restricted 
to the non-relativistic limit, i.e., it does not include the elfects of the 
existence of a marginally stable orbit of radius r ms . As mentioned by 
these authors, their model should be used only where ro 3> r ms . To es- 
tablish a more general model that can be applied to compact and mas- 
sive configurations near the mass-shedding limit, we have amended 
this shortcoming here by introducing Eq. (1121 . We demonstrate in 
this Appendix however, that in the case of moderately fast-spinning 
PSR J1903+0327 this refined approach is not essential. 

The upper panel of Fig. UOI shows three characteristic radii of the 
problem, ro, stellar radius R, and r ms , calculated along a spin-up 
evolutionary track. For the final configuration, / = 465 Hz, M = 
1.67 Mq, and B = 2 X 10 8 G, where ro reaches its minimum value, 
ro/r ms ~ 2. We compare two results, one where the effect of General 
Relativity is mitigated by Eq. 11211 . and the original model of KR, 
/ ms = 1 in Eq. lT2t . The difference in ro values is ~ 200 m, i.e., 0.7%. 
Lower panel shows the effective radial potential V (r, l(r)) (l(r) is the 
particle specific ang ular momentum on the circular orbit of the radius 
r; see, e.g., Eq. 2 of lBejger et al.ll2010l ); the final ro is sufficiently far 
from the r ma , hence we conclude that in the particular case of PSR 
J1903+0327 the relativistic correction necessary to account for r ms 
can be omitted. 
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